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Abstract 

We consider a simple model of the evolution of the concentration of a tracer, 
subject to a background shear flow by a fluid with viscosity u <C 1 in an infinite 
channel. Taylor observed in the 1950’s that, in such a setting, the tracer diffuses 
at a rate proportional to 1/iy rather than the expected rate proportional to v. 
We provide a mathematical explanation for this enhanced diffusion using a 
combination of Fourier analysis and center manifold theory. More precisely, 
we show that, while the high modes of the concentration decay exponentially, 
the low modes decay algebraically, but at an enhanced rate. Moreover, the 
behavior of the low modes is governed by finite-dimensional dynamics on an 
appropriate center manifold, which corresponds exactly to diffusion by a fluid 
with viscosity proportional to 1/za 


1 Introduction 


Taylor diffusion (or Taylor dispersion) describes an enhanced diffusion resulting from 
the shear in the background flow. First studied by Taylor in the 1950 ’s hie] many 
further authors have proposed refinements or extensions of this theory [31 lU El E]- 
In the present paper we show how in a simplified model of Taylor diffusion we can 
use center manifold theory to simply and rigorously predict the long-time behavior 
of the concentration of the tracer particle to any desired degree of accuracy. We 
note that in terms of prior work on this problem our approach is closest to that of 
Mercer and Roberts, |5j, who also use a formal center manifold to approximate the 
Taylor dispersion problem. However, they construct their center-manifold in Fourier 
space, an approach which is difficult to make rigorous because there is no spectral 
gap between the center directions and the stable directions. By introducing scaling 
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variables we show that a spectral gap is created which allows us to rigorously apply 
existing center-manifold theorems to analyze the asymptotic behavior of the problem. 
We believe that a similar approach will also apply to the full Taylor diffusion problem 
and plan to consider that case in future work. 

We consider the simplest situation in which Taylor dispersion is expected to occur, 
namely a channel with uniform cross-section and a shearing flow: 

d t u = vAu — A{1 + x(y))dx u ) — co < x < oo , — tt < y < n (1) 

u = u(x,y,t), u y (x, ±1, t) — 0. 

We assume that A is constant and that the background shear flow has been normalized 
so that x(y)dy = 0. Thus, the mean velocity of the background flow is A and 
we can transform to a moving frame of reference x = x — At. In this new frame of 
reference (and dropping the tilde’s to avoid cluttering the notation), we have 

d t u = v A u — Axd x u. 

We begin with a formal calculation that will be justified in an appropriate sense in 
subsequent sections and that provides some intuition about the expected behavior 
of (jl|. Given the geometry of this situation it makes sense expand u in terms of its 
y-Fourier series. Therefore, we write 

u(x,y,t) = '}Tu n {x,t)e my , and x(v) = ^2x n e tnv , 

n n 

where 

u n {x,t) = ^J e- %ny u(x,y,t)dy, x n = J e~ my x(y)dy, 

and we find (considering separately the case n — 0 and n^O) 

n = 0; d t u 0 = vd 2 x ii 0 - A(xu x ) 0 (2) 

n ± 0; d t u n = u(d 2 x - n 2 )u n - A(xu x ) n , 
where ____ __ 

( X u x) n = 'y ^ Xm(hn-m)x- 

m 

We now introduce scaling variables. These variables have often been used to analyze 
the asymptotic behavior of parabolic partial differential equations and they have the 
additional advantage that they frequently make it possible to apply invariant manifold 
theorems to these problems [7j. We expect that the modes with n / 0 will decay 
faster than those with n = 0, so we make a different scaling - of course we have to 
verify that the behavior of the solutions is consistent with this scaling. Let 

“ o(M) = 7m w ° {vTTt’ log(1 +1] 

Un(x,t ) = l +t ^ Wn + t )) ’ n ^°- ( 4 ) 
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In section )]3j below, we show that u n (x,t ) is basically an x derivative of a Gaussian, 
which generates the extra \ft decay. Proceeding, note that if we consider the advection 
term, the contribution to the n — 0 equation is of the form 


^ ^ Xm('d—m)x 
ra/0 




where we have no contribution from the term with m — 0 since \ m = 0 because x 
has zero average. The scaling in ([3]) and Q was chosen so that the terms in this sum 
will have the same prefactor in t as all the remaining terms in the equation for w o. 
More precisely, consider the various terms in the equation for uq. We find 


d t u 0 


11 11 
2(l+t)3/2 U ’° _ 2(1 + t)3/2 


£d£w 0 + 


1 

(1 + t)3/2 


d T w 0 , 


where we have introduced the new independent variables 

£ = 7 = , t = log(l + t) . 

VT+T 

Likewise, we have 

= 7TT7W^°’ ^ = 7T ,\ 3 / 2 ^ 

' ' m yo ' ' rr+0 

Thus, the equation for w 0 becomes 

d T w 0 = Cw 0 ~A^ y m (w_ m )g , (5) 

r7+ 0 


where 


£u> = z/(9|u’ + -d^w) . 

Remark 1.1. The spectrum of C can be explicitly computed. See tf^for more details. 


Repeating the calculation above for the evolution of the terms u n with n ^ 0, one 
finds that the terms are not any longer of the same order in t. Consider first the 
advective term which now has the form 

^ ' J Xm{'hn—m)x = ^ ^Xn^W o + + 3 / 2 ^ ; Xmd^W n - m . 

m ' ' m^n 

Working out the form of the remaining terms in ([2]), one finds 
( vr?w n + Axnd^wo) = e _r ((£ + 1/2)w n - d T w n ) - Ae~ r/2 ( ^ x m d^w n - m ) . (6) 

\m^n / 


The terms on the right hand side of this expression should go to zero exponentially 
fast so that in the limit r —>■ oo, w n satisfies the simple algebraic equation 

vn 2 w n + Axnd^wo = 0 . (7) 
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Remark 1.2. This is reminiscent of various geometric singular perturbation argu¬ 
ments and we will expand more upon this point in section 2. 


If we are interested in the long time behavior of the system, we can conclude from 
([tJ) that 

AXn 

-ydytco . 


W r, = 


V1W 


If we now insert this expression into (|5]) we find that 


A 2 \ 1 

d T w 0 = Tw 0 H-V -^XmX-md’tW o . 

v m z 5 

m^O 


Since x is real, X-m — Xm, and so the last term in the preceding equation can be 
rewritten as 

v m 2 I 

m^O / 

which implies that the equation for Wq becomes 


<9f w 0 , 


d T w 0 = (v + dfwo + ^<%(£w 0 ) • 


This is just the heat equation written in terms of scaling variables, but we see that 
the diffusion constant v has been replaced by the new diffusion constant v + Dt/v 
where the Taylor correction is 


d t = a 2 


\Xn 


rrt^O 


77r 


m 


That is to say, if we “undo” the change of variables, ([3]), and rewrite this equation i 
terms of the original variable uo(x,t), we find 

d t Uo = {u + D t / u)d 2 u 0 . 

Thus, we see that u 0 (which gives the average, cross-channel concentration of the 
tracer particle) evolves diffusively, but with a greatly enhanced diffusion coefficient. 

Remark 1.3. One can also check that the Taylor correction Dt to the diffusion rate 
computed above is the same as that given by the more traditional approaches cited 
earlier. 


In order to justify the above formal calculation, we need to analyze the system (|5]) 
<§ , which we rewrite here: 

d T w 0 = jCw 0 ~ A^2 Xm(w- m )z 


ra/0 


d T U)n 


(T + 1/2 )w n - Ae T/2 ( ^ Xmfyv-T-m J - e T ( vn 2 w n + Ax n d^w 0 ) 

\m^n / 
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In particular, we would need to show that there is a center manifold given approxi¬ 
mately by {w n = -(Axn/(vn 2 ))dzw 0 }. 

Rather than studying this full model of Taylor diffusion, we focus in this paper on 
the simplified model 

d T w = Cw — d^v 

d T v = (£ + l/2)v — e T (uv + d^w), (8) 

which corresponds to just the modes w 0 and w i (or more generally, to w 0 and u> n , 
where n is the first integer for which 7 ^ 0). The reader should note that this is 
not meant to be a physical model, but instead it is meant to be an analysis problem 
which reflects the core mathematical difficulties of analyzing the full Taylor Dispersion 
problem. Proceeding, note that the term proportional to e T//2 has disappeared since, 
if only wo and w\ are non-zero, that sum reduces to Xod^w \, and xo = 0. (We have 
also rescaled the variables so that the coefficient Ayi = Ax~ 1 = 1.) Also note that 
Q, written back in terms of the original variables, which we denote by w(x,t) and 
v(x,t), is given by 

W t = VW XX ~ V x 

v t = VVxx - vv - w x . (9) 

The classical picture of a center manifold would imply (see Figure [Tj) that solutions 
exponentially approach the invariant manifold (say at some rate a), with the dynam¬ 
ics on the center manifold given by the heat equation with the new Taylor diffusion 
coefficient. However, our analysis of system 0 . below, will show that the Taylor dif¬ 
fusion is really only affecting the lowest Fourier modes for w. Thus, the high Fourier 
modes still decay like e~ u ^ t for all \k\ > |&o|. Although this rate can be made uni¬ 
form for \k\ sufficiently large, it does not reflect the large diffusion coefficient of order 
1/u. If there really was a center manifold as suggested by the formal calculations, 
with dynamics on the manifold given by w t = (z/ + l/i/)w xx , then on that manifold all 
Fourier modes would decay like t . Note that this does not contradict the fact 

that Taylor diffusion seems to be observable in numerical and physical experiments 
(see, for example, the original experiments by Taylor IB. 0 ). The reason is that, for 
high modes with \k\ > (, the uniform exponential decay in Fourier space implies 

exponential in time decay in physical space, whereas low modes exhibit only algebraic 
temporal decay. Since algebraic decay, even relative to a large diffusion coefficient, 
is slower than exponential decay, the fact that the high modes do not experience 
Taylor dispersion does not prevent the overall decay from being enhanced by this 
phenomenon. 

While we believe that the picture sketched above applies to the full Taylor diffusion 
problem, as previously mentioned in this paper we analyze instead (J8]). For this 
coupled system of two partial differential equations we will show that 

• The long-time behavior of solutions can be computed to any degree of accuracy 
by the solution on a (finite-dimensional) invariant manifold. 
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• To leading order, the long-time behavior on this invariant manifold agrees with 
that given by a diffusion equation with the enhanced Taylor diffusion constant. 

• The expressions for the invariant manifolds can be computed quite explicitly, 
but we are not able to show that these expressions converge as the dimension 
of the manifold goes to infinity. Indeed, we believe on the basis of the argu¬ 
ment above, that (|8| probably does not have an infinite dimensional invariant 
manifold. 



Figure 1: Illustration of the invariant manifold one would expect for system , based 
upon the formal calculations. 


Our analysis will proceed as follows. First, in j|2j we’ll analyze the dynamics on 
the center manifold to see how the coupling between w and v affects the enhanced 
diffusion associated with the Taylor dispersion phenomenon. Intuitively, this is the 
key point of our result, and the latter sections can be thought of as justification of 
this calculation. Next, in j|3j we will obtain some a priori estimates on the solutions 
w and v in Fourier space and show that, to leading order, the long-time dynamics are 
determined by the low Fourier modes. Finally, in jj4| we show that the dynamics of 
these low Fourier modes are governed only by the dynamics on the center manifold 
analyzed in §[2j and thus, to leading order, the solutions exhibit the above-described 
enhanced diffusion. 

Our main result can be summarized in the following theorem, which is an abbreviation 
of Theorem 14.11 

Theorem 1.4. Given any M > 0, there exist integers m,N > 0 such that, for initial 
data {wq,vq) e L 2 (m) (see Definition \2.3\) , there exists a (2IV + 3) dimensional system 
of ordinary differential equations possessing an (N + 2) dimensional center manifold, 
such that the long-time asymptotics of solutions of ©, up to terms of Oft M ), is 
given by the restriction of solutions of this system of ODEs to its center manifold. 
Moreover, the dynamics on this center manifold correspond to enhanced diffusion 
proportional to v + \jv. 
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Remark 1.5. In the more precise version of Theorem TJ_ stated in Section 4, there 
-dependent constants that appear in the error term that make it clear that in 


are v- 


order to actually see Taylor Dispersion in the system, one has to wait at least a time 

t>0 (^A). 
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2 Dynamics on the Center Manifold 


In this section we focus on a center-manifold analysis of the model equation (| 8 ]). Our 
analysis justifies the formal lowest order approximation vv + d x w = 0 and shows that 
to this order the solutions behave as if w was the solution of a diffusion equation 
with “enhanced” diffusion coefficient Vt — (V+ p). Furthermore, the center-manifold 
machinery allows one to systematically (and rigorously) compute corrections to these 
leading order asymptotics to any order in time. 

Remark 2.1. As noted in the introduction, we do not expect that that the model 
equation (J 8 ]) (or the full Taylor dispersion equation) has an exact, infinite dimensional 
center-manifold. What we will actually prove is that, up to any inverse power of 
time, 0(t~ M ), there is a finite dimensional system of ordinary differential equations 
that approximates the solution of the PDE (J 8 ]) up to corrections of 0(t~ M ) and that 
this finite dimensional systems of ODE’s has a center-manifold with the properties 
described above. 

Because we expect v ~ - i.e. because we expect v to behave at least asymp¬ 

totically as a derivative, we define a new dependent variable u as 

v = d^u . ( 10 ) 

Inserting into the d T v equation in (| 8 ]), we get 

d T (d^u) = d T v — (C + 1/2) v — e T (vv + dguf) 

= (C + 1 / 2 ) d^u — e T (vd(U + d^w) 

= d^Cu — e T {yd^u + dyv) 

where we have used the fact that d^Cu = Cd^u + \d(U. After antidifferentiating the 
last line with respect to £, we get a system in terms of w and u: 

d T w = jCw — d^ 2 u 

d T u = Cu — e T (uu + w). ( 11 ) 
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Remark 2.2. Note that, if u G L 2 (m), the change of variables (10) implies that 
./- v(f, t)df = 0. We believe that, via minor modifications, our results can be ex¬ 
tended to the case when ff° v(f,t)df 0. We plan to discuss such modifications in 
a future work, when we study the full model g - g. 


Studies of Taylor dispersion generally focus on localized tracer distributions. For that 
reason, and also because of the spectral properties of the operators C which we discuss 
further below, it is convenient to work in weighted Hilbert spaces. 

Definition 2.3. The Hilbert space L 2 (m ) is defined as 

L 2 (m) = {/ € L 2 ( R) I ||/E = J (1 +er\m\ 2 d( < . 


Note that we require the solutions of the equation to lie in these weighted Hilbert 
spaces when expressed in terms of the scaling variables. If we revert to the original 
variables then it is appropriate to study them in the time-dependent norms obtained 
from these as follows: 



Thus, when we study solutions of our model equations in the “original” variables, as 
opposed to the scaling variables, we will also consider the weighted L 2 norms of the 
functions, but the different powers of x will be weighted by a corresponding (inverse) 
power of t to account for the relationship between space and time encapsulated in the 
definition of the scaling variables. These norms are discussed further in Section [3j 

Since we expect uu + w ~ 0, we further rewrite 0 by adding and subtracting ^<9|u> 
from the first equation and from the second finally obtaining 

d T w = C T w — - ( d^ 2 w + ud^ 2 u) 

d T u = Ctu - d^u — e T (uu + w ), (12) 

where 

C T f = (u + ^ + ^(£0) . 
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Thus, Ct is just the diffusion operator, written in terms of scaling variables, but with 
the enhanced, Taylor diffusion rate, ut — v + 1/v. 

The operators Ct have been analyzed in [9j. In particular, their spectrum can be 
computed in the weighted Hilbert spaces L 2 fm) and one finds 

a(C T ) = |A 6 C | SR(A) < j |} U {-t | k e n| . 


Furthermore, the eigenfunctions corresponding to the isolated eigenvalues X k = —/c/2 
are given by the Hermite functions 

MO = ^ / _L_ e~ g2/(4iyT) , and <j> k (g) = dgMO 

and the corresponding spectral projections are given by the Hermite polynomials 

H k (£) = 2 ^ e e/ ^ VT) d\e- e/ ^ VT) . 

Remark 2.4. The expressions in m for (f>k and H k are derived in the case when the 
diffusion coefficient is 1. The expressions given here follow easily by the change of 
variables f —* £/ y/uf. More explicitly, for the classical Hermite functions <f>o(y) = 
^_j=e -y2 / 4 , f k [y) = dyf o, and H k (y ) = |^e y 2 ' / 4 c^e -y2 / 4 ? one has the orthonormality 

relations f H k (y)4>e(y)dy = Sk/- Changing variables to y = tf/yfuf leads to the for¬ 
mulas for the eigenfunctions and spectral projections for Ct- Note further that with 
this definition, the Hilbert space adjoint of Ct satisfies Ct^H k = —| H k . 

Given the spectrum of Ct discussed above, we expect that the leading order part of 
the solution as t tends to infinity will be associated with the eigenspace corresponding 
to eigenvalues closest to zero. With this in mind, fix an integer N and assume that 
m > N + 1/2. This insures that the spectrum of Ct has at least IV + 1 isolated 
eigenvalues on the Hilbert space L 2 {m) and that the essential spectrum lies strictly 
to the left of the half-plane {A G C | 9?(A) < —IV/2}. Now define P/v to be the spectral 
projection onto the first N + 1 eigenmodes 

N 

P N w = 

k=0 


where 

a k ( T ) = 


We will write the solutions of (12) as 


{H k ,w(r)) L 2. 


w = Pnw + w s 
u = Pnu + u s . 


(13) 
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Based on the spectral picture and our discussion above, we expect that w s and u s 
will decay faster than Pnuj and Pnu (a fact which we demonstrate in Section [4]) and 
hence, since we are interested in the leading order terms in the long time behavior, 
we focus our attention on P^w and P^u. 


We will show that for any N the equations for P^w and Pn'u have an attractive center 
manifold and that the motion on this manifold reproduces and refines the expected 
Taylor diffusion. 


If we apply the projection operator P N to both of the equations in (12), we obtain 

N N k 1 N ~ 2 

dk^k = ^ — -otk&k ~ ~ ^ («fc + vf3k)<t>k+ 2 

fc=0 k=1 k=0 


N 


N 


N -2 


N 


y Mk = y - ^ Pkfa - - y Mk+ 2 - e r ( y + a k ) 4> k 

k=0 k=0 k=0 


^ k=0 


Shifting indices and matching coefficients gives us the following system of ODEs for 



and (d k '■ 

do 

= 0 


1 

b i 

~2 ai 


k 

dfc 

= -2 a* - 


-a k -2 + Pk -2 for 2 < k < N 


v 

$0 = ~e T (u/3q + « 0 ) 

= -^/3i - e T (vfii + ai) 

$k = ~ 7 }Pk - ~Pk -2 - e T {yf3 k + a k ) for 2 < k < N 


(14) 


Note that these equations contain no contributions from the “stable” modes w s and 
u s . Note further that, because of the form of the equations, those with even indices 
k decouple from those with k odd. Thus, we can analyze these two cases separately. 
We’ll provide the details for the case of k even below - the equations with k odd 
behave in a very similar fashion. 


Remark 2.5. Note that if we multiply all of the equations in (14) by e~ T and set 
e~ r = e (since we are interested in large times) we get equations that are formally of 
classical singularly perturbed form. (However, the small parameter e is time dependent 
here.) Invariant manifold theory has been a powerful tool in the rigorous analysis 
of singularly perturbed problems and that analogy will guide our use of the center- 
manifold theory in what follows. 


In order to make the invariant manifold more apparent we rewrite the even index 
equations by rescaling the time variable as 

r = log(l + t). (15) 
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In analogy with the above remark about singularly perturbed systems, we are es¬ 
sentially switching to a “fast” version of our system by making this change of time 
variable. Continuing, we introduce a new dependent variable 


rj — e T = 


1 

1 +1 


(16) 


Then, if we denote by a prime ', we have 

ckq = 0 

, (k 1 0 \ 

&k ~ ~ r l ( 2 ak "7 + P k ~ 2 ) 

Po — ~ Mo + «o) (17) 

Pk = ~ ( u Pk + «fe) - V Q/3fc + 2 J 

n' = -v 2 , 

where the values 2 < k < N are even. Notice the linearization of this system at the 
fixed point a k = p k = rj = 0 has eigenvalues A c = 0, with an [iV/2] + 2 dimensional 
eigenspace and A s = —u, with an [N/2] + 1 dimensional eigenspace (here [M\ refers 
to the greatest integer less than or equal to M). We proceed by diagonalizing the 
linear part of the system via 




&k 

—Olk + Pk 

V 


which transforms 0 into 


= 0 


Clu — 


b' = 


k 

-V ( g + bk ~ 2 


-ub 0 

-ub k 


V, 


' k 1 2 

V ( 77 bk -+ ~b k -2 


v* 


where again 2 < k < N are even. 


(18) 


(19) 


The remainder of this section is devoted to the analysis of these equations and we 
prove two main results: 


We first show that, for any N, (19) has a center-manifold of the type described 
in the introduction, and we derive explicit expressions for the functions whose 


graphs give the manifold. (See Propositions 2.6 and 2.7 
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We derive the asymptotic (in r) behavior of solutions of these equations. (See 


Propositions 2.8 and 2.10 and Corollary 2.9 


We begin by noting that the linearization of (19) at the fixed point a*, = bk — r) — 0 
has eigenvalues A c = 0, with an [1V/2] + 2 dimensional eigenspace and A s = — z/, 
with an [1V/2] + 1 dimensional eigenspace. Thus, from the classical center-manifold 
theorem (say, for example, the center-manifold theorem proven in [[TOj), we know 
that (at least in a neighborhood of this point), there will be an invariant [N/ 2] + 2 
dimensional center manifold. We also know that, in a neighborhood of the origin, the 
center-manifold can be written as the graph of a function with components 


bk — hk{ciN, ..., a 0 , i])- 


( 20 ) 


In addition, because of the “lower triangular” form of the equations (i.e. the fact that 
the equations for a' k and b' k depend only on ae and be with £ < k), we find that we 
can express the manifold as 


bk — h k (a k , ak- 2 ,... ,ao,rj) . 

We now show that we can find explicit expressions for the functions h}. successively, 
starting with ho and then progressing through h 2 , h 4 , etc. What’s more, these ex¬ 
pressions hold for all ak,CLk- 2 ,... ,ao,rj, i.e. without the restriction to a small neigh¬ 
borhood that is inherent in general center-manifold theorems like that of [ 10 ]. 

We start with the equations for ao and b 0 which are just 

a' 0 = 0 
b' 0 = -vbo . 

From this we see immediately that we can choose the invariant manifold to be the 
graph of ho = 0. However, note that this example also reminds us that the center 
manifold is not unique, since we could also choose the center manifold to be given 
by the graph of ho(«o ; 'h) = A’oe _ 1 '' /j? ao- This is consistent with the theorems on the 
existence of center manifolds, since both of these manifolds have the same Taylor 
expansion to any finite order about ao = V = 0. For simplicity, in what follows we 
will always use the first function - i.e. we will take ho = 0. 

Now consider the center manifold for a 2 and b 2 . Since the equations for a 0 , a 2 , b 0 , b 2 , r] 
decouple from all other a& and b k , we expect the center manifold to be given by the 
graph of a function b 2 = h 2 (a 2 , a 0 , ?/). In fact, as we show below, it has no dependence 
on a 2 - i.e. we can take b 2 = h 2 (ao, ?/). In this case the equation for the invariance of 
the graph of this function takes the form 

11 2 r/ 1 

(D ao h 2 )a'o + (D v h 2 )rj' = -rjh 2 - uh 2 - —ho + —^h a o ■ 
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Inserting the equations for a' Q and rj' and using the fact that ho = 0, we find 

X 

-r] 2 (D v h 2 ) = — 77/12 - vh 2 + r/ao ■ 


We now show that h 2 is linear in Oq, so we write 


h 2 (oo, h) = 02,0(h)«o , 

and find 

-hV 2 ,o = ~(V + *002,0 + • 

This equation is hard to solve in general due to the singular point at 77 = 0, but 
remarkably, 

0 2 , 0 ( 77 ) = -3 

is an exact solution (which goes to zero as rj —* 0 ), so 

h 2{ao,V) — —y 

is a function whose graph (together with that of h 0 = 0 ) gives us the center manifold 
for the equations for ao, a 2 , b 0 , b 2 , 77 . Due to the singular point at 77 = 0, this may not 
be the only solution (just as in the case for ho), but we are free to choose this special 
solution for h 2 . 


Next we consider the case of ^, 4 ( 04 , a 2 , ao, rj). Building on the examples above we 
show that 


• h 4 is independent of a 4 ; 

• h 4 is linear in a 2 and a 0 . 

If this is the case we can write 

h^{a 2 , Oo, T/) = 04,2(h)®2 + 04,0(h) a o • 

Inserting this form of the solution into the equation for the center-manifold, we find 

04,2 faK + 04,0 (h) a 0 + (04,2 (h)®2 + <f>4fl(v)aoW 

= -{y + 277 ) (04,2 ( 77)02 + 04,o(h) a o) + - — 1 -^r 

where in the last term we have plugged in the expression for h 2 . Inserting the equa¬ 
tions for o!. 2 and 77 ' and grouping the terms proportional to a 2 and ao we find two 
ODE’s for the 0's, namely 

-h 2 04 , 2 (h) = -iy + 77)04,2(77) + 4 

2 / 7 2 

-h 2 04,o(h) = -O' + 2/7)04,0(77) - . 
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The first of these equations is the same as the equation for 02,o above so we have 

04,2 (V) = ^ • 

The second equation is very similar and we find that it again has a simple, exact 
solution, namely 

i ( \ 2ry 

Mv) = - /y5 • 

Tims, we also have an exact expression for the center-manifold in this case: 

, / x VO -2 2 T] 2 a 0 

n 4 (a 2 , a 0 ,//) — —~- . 


One can continue this procedure. For instance, for the function h 6 , 
formula 


hs{a 4 , a 2 , ao; v) 


rja 4 




2 ifa -2 5r/ 3 a 0 

v 5 v 7 


one obtains the 


This leads to the following 

Proposition 2.6. For any k — 0, 2,4,..., there exist constants {H(k, k — 2£)} such 
that the graph of the function 

k/2 

h k (a k - 2 , a k - 4 , ■ ■ ■, a 0 , rj) = y H(k, k - 2k)rf a k _ n (21) 

t=\ 

gives the invariant manifold for b k . Furthermore, for any fixed k, the coefficients 
{H(k, k—2£)} can be explicitly determined, and the coefficients H(k,p ) 


Proof. The proof proceeds inductively. Note that we have already verified the induc¬ 
tive hypothesis for k = 0,2,4. (We take the empty sum that occurs on the RHS of 
(21) when k = 0 to correspond to ho = 0.) Assume that it holds for all even integers 
less than or equal to k — 2. We now show that it holds for h k . 


Inserting our inductive hypothesis into the invariance equation we find 

k/2 k/2 

H(k , k - 2£)ri i a' k _ 2£ + ^ £H(k, k - 2F)rf~ 1 a k ^ 2 (X\ 

i =i i=i 

k 2 1 

= - 7 Pl h k - vh k - r)h k _ 2 + -^r]a k _ 2 


k/2 


k 


V 

A+i 


v* 


k/2 


y -H(k , k - 2£)rj l+l a k - 2 i - y vH{k, k - 2t)rf 


a k -2C 


i =\ 


e=i 


k/2-1 


— y H(k - 2, k - 2 - 2 t)rj e+1 a k - 2 - 2 t + —^V^k- 2 


£=1 


( 22 ) 
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Inserting the equations for a' k _ 2e and rf into the first line of (22), one finds 

k/2 


k/2 

Y H{k, k - 2 t)rf 

t=\ 

k/2 k 


k-2£ 


V a k-2£ — Vhk-2£-2 

k/2 


Y^{k,k-2£)rf 


£+1 


1= 1 


- Y ~H{k, k - 2l)rf+ x a k _ n - Y H(k, k - 2£)rf +1 h. 


k-2t-2- 


&k-2t 


(23) 


1=1 


1=1 


Note that the first sum in the last line of (23) cancels the first sum on the RHS of 
(22). Thus, we can rewrite (22)-(23) as 


k/2 

Y u H(k, k - 2£)p £ a k _ 2e 
e=i 


1 2 fc/2_1 . 

~o Va k -2 - Y -2,k-2£- 2)rj e+1 a k _ 2 -2£ 

u z v z —' 

1=1 

k/2 

+ Y k - 2£)ii +l h k -n- 2 . (24) 

i=i 


We now rewrite the last sum in this expression by using the inductive form of h k - 2 e- 2 , 

k/2-(£+l) 

h k ~ 2 £~ 2 = Y^ H{k — 2{i + 1), k — 2{l + m + l))p m a k - 2 ^ +m+ \) . 

m =1 


Thus, 


k/2 

YH(k,k-2£) V £+1 h k _ 2e _2 

t=\ 

k/2 k/ 2-(£+l) 

= E Z V(k, k - 2e)H(k - 2{i + 1), k ~ 2{£ + m + l)W +m+1 a k -2(£ +m+ i) 

1=1 m= 1 

k/2 p—2 

= EE H(k, k - 2£)H(k -2{£+l),k- 2p)p p a k . 2p , 

p =3 £= 1 


where in the last term we set p — i + m +1 and interchanged the order of summation. 
If in the last sum in the first line of (24) we also change the summation variable to 
p = £ + 1 we find that (24) can finally be rewritten as 


k/2 12 fc/2 

Y vH^k, k - 2£)ifa k - 2 i = —^V a k -2 ~ ~Y ~ ~ 2 p))rfa k ~ 2p 

1=1 U p= 2 

k/2 p—2 

+EE H(k, k - 2£)H(k -2(£+l),k- 2p)r/ p a k _ 2p . (25) 

p =3 £=l 
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We solve (25) for H(k, k — 2t) 1 beginning with H(k, k — 2 ). Since the only term on the 
RHS of (25) proportional to a *-2 is the first term, and we obtain H(k,k — 2) = dj, 
consistent with the inductive hypothesis. Next consider H(k,k — 4). In this case, 
we consider all terms in (25) proportional to Gq._ 4 . The only one comes from the 
second term on the RHS of the equation and we have H(k,k — 4) = — ^ H{k — 
2, k — 4). The inductive hypothesis implies that H(k — 2,k — 4) ~ so we 

find H(k, k — 4) ~ 0(v~ 5 ) as required by the inductive hypothesis. We now continue 
to solve for the coefficients H(k, k — 20), t — 3,4,..., noting that in each case, the 
terms on the RHS of the equation proportional to ak- 2 e have coefficients that have 
already been determined at prior stages of the inductive process and that they are all 
O^- 21 - 1 ) = 3 ). □ 


We now describe the entirely analogous results for the modes ak and (3k with k odd. 
If we introduce new variables t and tj as in (15), (16), and diagonalize the linear part 
of the resulting equations using the change of variables (18), we find: 


a 'k 

K 

K 


~ \^ 


-v 



a k + b 



v + ~V ) h 


-vb k - rj 



1 

2 + 
v l 




1 


where the values 3 < k < N are odd this time. 


(26) 


Proceeding as before, consider first the equations for cq, b±, and 7] which decouple 
from all the rest of the equations. Then by inspection we see that, just as for 6 0) the 
graph of the function h\(a\,rj) = 0 is an invariant center manifold for these equations. 
We now include the equations for 03 and 63 and, building on the experience from the 
even case, look for an invariant manifold of the form 


b 3 = h 3 (a 1 , i q) = 03,i ( 77)01 . 

Inserting this into the equations, we see that in order for this graph to be invariant, 
03 ,i must satisfy 

3 Tj 2 JJ 

03,1 a 1 + ai03,i rf = -(v + — 77)0 3 ,iai + - -jh x . 

From the fact that h\ = 0 and the equation for a \, we see that this reduces to the 
ODE for 03.1 

-^ 2 03 ,i = -iy + ^7)03,1 + 4 • 
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This is the same equation satisfied by 02,o and thus we find 

h3{ai,rj) — — t " . 

v 6 

Proceeding now as in the even case, we establish the following proposition by induc¬ 
tion. 

Proposition 2.7. For any k = 1,3,5,..., there exist constants {H odd (k, k — 2£)} 
such that the graph of the function 


k -1 
2 

h k (a k - 2 , cifc- 4 ,... , 01 , 77 ) = H odd (k, k - 2£)rfa k _ n 

1 =1 

gives the equation for the invariant manifold for b k . Furthermore, for any fixed k, 
the coefficients {H odd (k,k — 2£)} can be explicitly determined, and the coefficients 
H odd (k,p) ~ O^v-^-pI- 1 ). 


We conclude this section by using our expressions for the center-manifold to derive 
the asymptotic behavior of the coefficient functions a k and b k (or equivalently a k and 

M 

Begin by noting that from the general theory of center-manifolds, any solution with 
initial conditions in a neighborhood of the invariant manifold will approach the man¬ 
ifold at a rate 0(e- ut ) = 0{e~^ eT -^). Tims, we can determine the long time 
asymptotics of all solutions in this neighborhood by focusing on the behavior of solu¬ 
tions on the invariant manifold. Note that this means, for solutions with sufficiently 
small initial conditions, that after a time r such that ve T 1 , we will be very close to 
the center-manifold and the behavior of solutions on this manifold will determine the 
asymptotic behavior of solutions after this time. Reverting from our rescaled time r 
to the original time t in the problem this means that solutions on the center-manifold 
will determine the behavior of solutions for times t > O(^), which is the expected 
timescale for Taylor Dispersion to occur. At the moment, it appears our results only 
hold for solutions with small initial conditions. However, it turns out our formulas 
for the center manifolds (which are defined globally) are also globally attracting on 
the timescale t > O(^). We provide details in the Appendix. 

We proceed with our calculation of the asymptotics of the quantities a k and b k . As 
in the case of the calculation of the manifold we focus separately on the coefficients 
with even and odd indices. Starting with the coefficients with k even, note that we 
obviously have «o = constant, so we begin with k — 2. 


Given 


a 2 = - 77(02 + b 0 ), 


we can simplify this by noting that bo = ho = 0 on the center-manifold. Finally, 
it’s simpler to solve this differential equation by reverting from the t variables to 


r = log(l + i); keeping Remark 2.5 about singularly perturbed systems in mind, 
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notice we are essentially switching to the “slow” version of the system (which gives 
the dynamics on the center manifold). The equation then reduces to 

CL2 = —CL2 , 

from which we can immediately conclude that 

a 2 (r) ~ 0(e~ T ) . 

Next consider a^, for which we have (again, rewriting things in terms of the temporal 
variable r) 

o u o e ~ Ta ° 

CL 4 — —ZO4 — 0 2 — —ZCI4 — --— , 

V 6 

where the last equality used the fact that b 2 = h 2 (ao, rj) — ^ on the center-manifold. 
Solving this equation using the method of variation of constants, we find that 

<* 4 (r) ~ Of^) ■ 


As a last explicit example, consider the case of a6 where we have 
clq — ^4 — 3 clq 


e r a 2 2e 2r a 0 


Finally, since a 0 is constant and a 2 (r) ~ 0(e T ), we see that the asymptotic behavior 
of a 6 is 

06(r) ~ of^r-) • 

We can generalize these results in the following 


Proposition 2.8. Suppose k — 4, 6 ,... is an even, positive integer. On the center 
manifold of the system of equations (19), the variables a*, have the following asymp¬ 
totic behavior: 


a k (r) | < 


C(N,k)e~% T 

jyk—l 

fc+2^ 

C(N,k)e —3~ r 


l^k — l 


: k = 0 mod 4 
: k = 2 mod 4 


Note that once we have these formulas, the expressions for the center-manifold im¬ 
mediately imply the following. 


Corollary 2.9. Suppose k — 4, 6 ,... is an even, positive integer. On the center man¬ 
ifold of the system of equations (19), the variables b k have the following asymptotic 
behavior: 


Mr) | < 


C(N '%^ T : k = 0 mod 4 

fc +2 

C(N %^ T : k = 2 mod 4 
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Proof. The proof of Proposition [278] is a straightforward induction argument. Suppose 
that we have demonstrated that the estimates hold for k — 4, 6 ,..., ko- We then show 
that it holds for ko + 2. The equation of motion for a,k 0 +2 is 


_ k 0 + 2^ u 

^fco+2 q ^/co+2 '^ko\Q j ko— 2? &ko— 4? • • • i ^ /• 


Inserting the formula for /i fco from Proposition 2.6 and solving using Duhamel’s for¬ 
mula, we obtain the bound 


K 0 + 2 I < 


C(N) 


fco/2 

E 

e=i 


a ko—2l n £' 


(27) 


Consider the case fco = 0 mod 4. Then 


f 2 mod 4 if £ is odd 
( 0 mod 4 if i is even 


and correspondingly from the induction hypothesis, 


Wk 0 -2e\ < 


C(AQe~ 


— fc 0 —2^+2 


y fc-i 

-kr,-2£ 

C(N)e~ 4 ^ 

j.k — l 


if i is odd 
if t is even. 


Inserting into (27), using the fact that 77 = e T , and splitting the sum into even and 
odd we obtain 


^fco +2 


C(N) 

v 


fco/ 2-1 

E 


1=1 ,£odd 


_ (fc 0 -2C+2)r 

e 4 e 

jyko—2£—lj^2£ 


feo/2-2 (k 0 -2t)r 


It 


+ 


V kv-21-\ V 2£ 

1=2 ,teven 


k Q T 

a 0 e 2 


(28) 


Notice we have to separate out the i = k 0 /2 term because this corresponds to ao, 
which is actually constant. We are interested in locating the slowest decaying terms. 
These terms will have, in the exponent, the least negative coefficients on r. For £ > 1 
odd, the coefficients in the exponent are 


ko-2i + 2 ko 1 t 

4 ’ _ ~ 2~ 2 


(29) 


which are least negative when i = 1. The corresponding coefficient in the exponent 
is — and so the slowest decaying term from the i odd sum is 0(e~~^ T ). We 
determine the slowest decaying term in the t even sum. For i > 2 even, the coefficients 
in the exponent are 



(30) 


which are least negative when ^ = 2. The corresponding coefficient in the exponent is 
again — and so the slowest decaying term from the l odd sum is again 0(e~~^ T ). 
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Lastly, we determine the v dependence of the constant. The largest power of v in the 
denominator comes from I = fco/2 and is ^ +1 . Therefore we have 


a k 0 +2 


< 


C(N) _ko±± T 

—-—-e 4 

iyko~\~l 


Recalling that we are in the case k 0 = 0 mod 4 ( so that k 0 + 2 = 2 mod 4), we 
have verified the claim in this case. The case ko = 2 mod 4 follows similarly. Once 
Proposition |2.8| is established, a nearly identical calculation establishes Corollary 24) 


□ 


The coefficients a*, and bk, with k odd, can be estimated in an entirely analogous 
fashion to obtain the following proposition. 


Proposition 2.10. Suppose k = 1,3, 


manifold of the system of equations (26), 
totic behavior: 


is an odd, positive integer. On the center 
the variables ak have the following asymp- 


Ofc(r)| < 


C(N,k)e 


. fc+i. 


C(N,k)e 


. k+ 3, 


: k = 1 mod 4 
: k = 3 mod 4. 


(31) 


If k = 3, 5,... (recall that b\ = 0 on t/ie center manifold), the corresponding coeffi¬ 
cients b k satisfy the estimates 


\h(r) | < 


C(N,k)e~ 

isk+i 

C(N,k)e~ 

l^k+l 


5+fc 


3+fc 


: fc = 1 mod 4 
\ k = 3 mod 4. 


(32) 


3 A priori estimates via the Fourier Transform 


In order to show that the center manifold, discussed in the previous section, really 
does describe the leading order large-time behavior of solutions of (|9]), we need to 
make our discussion before Theorem 1.4 in the introduction more precise (which 
basically says Taylor Dispersion only happens for low wavenumbers). We’ll have to 
undo the scaling variables, and switch to the Fourier side; this way we can precisely 
cut-off wavenumbers larger than, say \k 0 \ | and quantify how fast these “high” 

wavenumber terms decay. To do this in a way that is consistent with the analysis in 
section 2, we need to introduce a new norm ||| • |||, which, when applied to functions 
on the Fourier side, is equivalent to the L 2 (m) norm applied to their real-space scaling 
variables counterparts. 


The main result (see Theorem 4.1 in Section 4) depends on estimates of the solution 
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in L 2 (m). With this in mind, we note that 


\w(t) ||l 2 m < C(m)(t + 1) 1/4 


\ 


E 

3=0 


(i + ty/z 


diw{-,t) 


=■ Ill-Will, 


L 2 


and below we will bound each partial derivative of w(k,t). Note that || ■ ||L 2 (m) and 
|| ■ HI are indeed equivalent norms, which follows from the fact that 


which in turn implies that |||w}(t)||| < C(m)||w(r)||£ 2 ( m ). 


lL 2 (i)> 


Consider equation ([9]). Let w = Tw and v = Tv , where T sends a function to its 
Fourier transform. We obtain 


The solution to this equation is 


A(k) = 


—vk 2 —i k 
—i k —v(k 2 + 1)1' 


w(k,t) \ _ A(k)t 

v(k,t ) 1 


W o(k) 

vo (k) 




W 0 (x) 

Vo(x) 


To understand these solutions, we must understand which we’ll do by diago¬ 

nalizing A(k). The eigenvalues of A are given by 

A ±(k, v) = -vk 2 -^ + 7 >Vv 2 - 4 k 2 , 


and the corresponding eigenvectors are 

i k 


v±(\,k) = 


-vk 2 — X±(k, v) 


i k 

^ + Wv 2 -4k 2 ) ' 


We put these into the columns of a matrix S = [u+,u_] and obtain 

S = 


i k i k 

\[v — y/v 2 — 4k 2 ] \[v + y/v 2 — 4 k 2 ] 


S~ = 


1 


\[v + y/v 2 — 4 k 2 ] 


-i k 


i ky/v 2 - Ak 2 V|[ _z/ + Vv 2 - 4k 2 ] i k 


We then have A = SAS , where A = diag(A + , A_). 

Remark 3.1. Note that S becomes singular when k = ±zz/2, because for that value of 
k there is a double eigenvalue, and a slightly different decomposition of A, reflecting 
the resultant Jordan block structure, is necessary. This will be dealt with in the proof of 
Propositio7i \3.3[ We do not highlight this issue in the below formulas for the solution, 
as we wish to focus on the intuition for how to decompose solutions, which does not 
depend on this singularity. 
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Hence, 


^A(k,u)i = ^ 


,\ + (k,u)t 

0 


0 

_ (k,is)t 


S~ l (k,u), 


or explicitly 


w(k,t ) = 


i/c(e A_t — e A+t ) „ If (—v + y/v 2 — 4 k 2 )e x ~ t + (u + y/v 2 — 4 k 2 )e x+t 


y/v 2 — Ak 2 ° 2 


y/v 2 — Ak 2 


w 0 


*>(M) = - 


1 / (—v + y/v 2 — A k 2 )e x+t + (v + y/v 2 — Ak 2 )e x ~ i 


ik(e 


A-(-t _ — t 


y/v 2 — Ak 2 

Wo, 


Vo 


y/v 2 — Ak 2 

which we’ll abbreviate as 

w(k, t ) = (fi(k)wo(k) + h(k)v 0 (k)) e A+<fc) * + 9 (i)e'- ( ‘ ) ‘ 


(33) 


(34) 


and similarly for v. The motivation for separating the solution in this way is the fact 
that Re(A_(fc)) < — v/2, and so any component of the solution that includes a factor 
of e A- ( fc )* will decay exponentially in time, even for k near zero. Hence, it is primarily 
the first term, above, involving e A+ ( fc ^ that we must focus our attention on. We’ll 
proceed with the analysis only for w; all of the results for v are analogous. 

Remark 3.2. In order to justify the difference of (t + l) -1 / 2 in the scaling variables 
for w and v, corresponding to (|3]), we need to show that v decays faster than w by 
this amount. This can be seen from the above expression for solutions. In particular, 
for k near zero, say \k\ < v/2, we have 


D A(k,u)t 


n —isk 2 t 


ikv 


' „2 


An extra factor of k corresponds to an x-derivative, and so the v component does 
decay faster by a factor of £ -1 / 2 . 


We will split the analysis into “high” and “low” frequencies using a cutoff function 
and Taylor expansion about k = 0. Define 

n< = {|fc| < 

and let ip{k) be a smooth cutoff function equalling 1 on Q, < and zero for \k\ > ^-j^+z/ 2 . 
We then write w as 

w(k,t ) = ip(k)w(k, t) + (1 — ip(k)) w(k, t) 

=: ip(k) (fi(k)w 0 (k) + f 2 (k)v 0 (k)) e x+{k)t + if{k)g{k)e x ~ {k)t + w high (k , t). 
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Again, the motivation is to focus on the part of the solution that does not decay 
exponentially in time. This does not necessarily occur if k is small, which is exactly 
where if(k) ^ 0. 

Notice that, on f2 < , we can write \+(k) = — (y + k 2 + A(k), where 


m = 



(35) 


In h2 < , 4 k 2 /v 2 < 15/16 < 1, and so the above series is convergent. It will also be 
important that it starts with four powers of k. More precisely, 


m = 



This representation for A (k) holds by similar reasoning whenever if(k) ^ 0. We now 
write 

w(k,t) = ^{k)e~ UTkH e K{ - k)t {h{k)w 0 {k) + f 2 (k)v 0 (k)) + ^(k)g(k)e x ~ (k)t + w high (k,t) 
=: 'if(k)e~ UTk2t w(k, t ) + k)g(k)e x ~ [k)t + w high (k, t), 

where Ur — v + - and 

J V 

w(k,t ) = e A{k)t (/i(k)wo(k) + f 2 (k)v 0 (k)). (36) 


The purpose of this last part of our decomposition of solutions is to emphasize that, 
to leading order, the decay of the low modes will be determined by the term e~ l " rk t . 
Therefore, the Taylor dispersion phenomenon is also apparent in Fourier space. 


Finally, we Taylor expand the quantity w into a polynomial of degree N, plus a 
remainder term: 


w 


= Mjt* + 


3=0 


N 


w(k,t) - Y 


d 3 k w(0,t) 


3=0 


f- 


k? 


w, 


N 

low 


+ W 


res 
low ’ 


Thus, we have (suppressing some of the k and t dependence for notational conve¬ 
nience) 


w(k, t) = + wfZ) + i> S e X - m ‘ + (37) 

The main results of this section are 

Proposition 3.3. There exists a constant C, independent of u and the initial data, 
such that 


\\& k Whigh\\v + ||f>/!(WF A ^)|l/> 1 < Cv 2 J e ^(Ht&ollci + IMb)- 
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Proposition 3.4. There exists a constant C such that 





< 


L 2 


jv , j N, 1 (ll^ollc^+J + ll^ollc^+J')- 

V 4 + 2t 4+2 


T/ie constant C depends on N, but it is independent of v. 


Here \\f\\cj = Si=o su PfceR \®kf(k)\- These results imply that Whigh and gfje^ decay 
exponentially in t, and are thus higher-order, while decays algebraically, at a 
rate that can be made large by choosing N (which will correspond to the dimension 
of the center manifold from Section [2]) large. In the next section, §[4| it will be shown 
that the behavior of the remaining term, wjf )W , is governed by the dynamics on the 
center manifold, in which one can directly observe the Taylor dispersion phenomenon. 

Proof of Proposition |3.3| 

Notice that, for k G fP (the support of Whigh ), the eigenvalues A ±{k) both lie in a 
sector with vertex at (ReA, ImA) = (—ok 2 — z//4, 0). Therefore, to obtain the desired 
bound, we need to determine the effect of the derivatives <9j[,. Such a derivative could 
potentially be problematic, due to the factors of \/v 2 — 4 k 2 , which can be zero in fP. 
(This is exactly due to the Jordan block structure at k — ±u/2.) To work around 
this, we use the fact that we can equivalently write 

(w{k,t)\ _ A (k)t fw Q (k)\ 

\v(k,t)J \v 0 (k) J 

and bound derivatives of this expression for k G fP. Such derivatives either fall on 
the initial conditions, which leads to the dependence of the constant on the C norms 
of do and Po ; o r the derivatives can fall on the exponential. In the latter case, using 
the fact that 



which behaves no worse that 0(k), we obtain terms of the form (writing U = (wq, ho) 
for convenience) 

\\(kt) p e AWt d q k Uo \\ 2 L 2 <C\\U 0 \\ 2 Cg J \kt\ 2p \\e A{ - k)t \\ 2 dk. 

Next, note that ||e A ^ fe ^|| < Cv~ 2 e~ u ^ kZ+l / A " >t . This follows essentially from the above- 
mentioned bound on the real part of A± in fl 51 . One needs to be a bit careful when k = 
±z//2, as there A + = A_. This changes the bound from e A+ ( fc ) t to ~ ute A +( fc )i, but 
this power of t can be absorbed into the exponential since Re(A + ) < — uk 2 — v/A — dv 
for some 5 > 0 that is independent of u. The factor of v~ 2 that appears is related to 
the fact that ||S'~ 1 || = 0(v~ 2 ) for k G IP, k ^ ±u/2. Thus, we have 

\\{kt) p e A ( k)t d q k U 0 \\ 2 L2 < Cu- 4 \\U 0 \\ 2 Cq J \kt\ 2p e~ 2 ^ k2+u/4)t dk 

< Cv- A - p - 1/2 \\U 0 \\ 2 Cq t p - 1/2 e- ut/2 

< Cu- 4 - 2p e-^\\U 0 \\ 2 Cq , 
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which proves the result for Whigh■ A similar proof works for the ||c^('i/h?e A ~ f )||£2 term. 

□ 


Proof of Proposition |3.4| 


We now derive bounds on the residual term ife Urk2t w^. Recall the integral formula 
for the Taylor Remainder: 


«(M) 



rfcjv 


d k*^ m ( k n+^ t)dk„ +l dt„ ... 4 ,- 


(38) 


With this formula in mind, we want to derive bounds on the derivatives of w(k,t), 
but we need only deal with k G fU, since we are ultimately estimating the size of 


ipe 


—UTk 2 tyjres 


low’ 


Recall that 

w = e A{k)t ( fiw 0 + / 2 W) • 

The functions /1 and / 2 are smooth in 0 < , so our estimate will depend on derivatives 
of the initial data and derivatives of e A ^*. However, the reader should note that 
fi and / 2 are also dependent on z/, but their derivatives give us inverse powers of v 
no worse than any others appearing in this section, so we choose not to explicitly 
keep track of these powers. The following lemma will be used in estimating these 
derivatives: 

Lemma 3.5. Let = k d e~ UTk ~' t . Then 

\M;t)\\L*<C(d)(v T t)- 2 -^ 


Proof. Use the fact that f R e x2 / 4 dx = 2^/n and change variables. □ 

With this lemma in mind, we need to keep track of the powers of k and t that appear in 
d^ e Mk)t_ ^ see w j 1 y one expect the powers of v and t appearing in Proposition 

3.41 consider the following formal calculation. Recall from the Taylor expansion of 

_ - fc4 i 

A (k), we have w ~ e ^ . 


We are essentially estimating 


11^ 


3 — uxk-t —res 


W. 


low \\L 2 i 


with the aid of the estimate 

\\k d e - v T k2t \\ L 2 <C(d) (v T t)~*~* 


and the Taylor Remainder formula 


™ r iZ(k,t) 




JV+l) t)d kN+1 d kN . 


■ ■ d kl ■ 


(39) 
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1 — fe4 1 - 

We’ll proceed by finding bounds on a k e ^ , and plug into (39) with J = iV + 1. We’ll 

make the following changes of variable: we set 


x 


T = — 

i/ s 

= T 1/A k 


(40) 


so that 



and 

3 J k w = T J/A d J x w. 


(41) 


Let’s proceed by computing x —derivatives of w, only taking into account what powers 
of x appear at each stage. In the following, a prime means d x . We compute 

w' ~ x 3 e~ x4 
w" ~ ( x 2 + x 6 ) e~ x4 
w'" ~ [x + x 5 + x 9 ) e~ x4 . 


In particular, notice that the powers of x that appear in the J —th derivative can be 
obtained from the powers of x that appear in the J — 1st derivative by subtracting 
one from each power appearing (where only nonnegative powers are permitted), and 
also adding three to each power appearing: 

w (A) ~ (x° + x A + x 8 + x 12 ) e~ xi 
(x 3 + x 7 + x n + x 15 )e-* 4 
w^ ~ (x 2 + x 6 + x 10 + x 14 + x 18 ) e~ x4 . 


In general, we have 


J—2 


d * w ~ 


x R+Al e- x4 


1=0 


where R = (— J ) mod 4. In the original variables, we have, using (41) and (40), 




w 2 


R+Al 


e , 


or more precisely, 




R+il 


t 

e v 3 . 





Combining with the Taylor Remainder formula and setting J = N + 1, we have 




N - 1 Al/4)(R+N+l)+l 

< E ^3/4„K +1 V +1 ) + 3, ll^ +1 + B+4 ' e '" Tt s iu- 


2=0 


Using the estimate (3.5), we get 


N -1 


Me-" kH wlZ\\ L ’ < C(N)Y _; 


2=0 


l/4i?—1/4(7V+1)—2—1/4 
z/ l/4fl+l/4(JV+l)+Z-l/4 


= C(JV) ( , () - 1/4 ' i - 1/4< ' v+1) tU £ (iz-'r 1 ) 1 


2=0 


= C(N)(i/i) 


—1/4R—1/4(AT+1) 


t-1/4 / 1 — ' 


z/- 1 / 4 \ 1 - (zz-U- 1 ) 

Therefore if t > we have 

4- —1/4 

< 2C(JV) {vt )-V* R -mM _____ 

which implies that 




C(JV) 

JV JV I 1 
Z/ 4 t 4 _l_ 2 


as reflected in Proposition |3.4[ This concludes the formal calculation. We proceed 
with deriving the precise estimate. 

Because k is small in hU, powers of k are helpful, so we only need to record the 
smallest power of k relative to the largest power of t. We obtain additional powers 
of t when a derivative falls on the exponential (as opposed to any factors in front of 
it), which creates not only powers of t but powers of (A'(k)t). When derivatives fall 
on factors of A '(k) in front of the exponential, we obtain fewer powers of k but no 
additional powers of t. Using (35), we see that A '(k) ~ k 3 / z/ 3 , and so d J k e A< '- k ^ t will 
lead to terms of the form 


u 3 J y u 3 

This implies that 


umv /,yp V)1 


q + 2 li + 3^2 + 4/3 — j. 


(k 3 t\ q fk 2 t\ h fkt\ h f t 


V z/ 3 J \u 3 


v 3 ) U 3 1 6 


A (k)t 


\diw(k,t)\ < C(\\w 0 \\ C i + ||^oHcu 

for any q + 2li + 3Z 2 + 4Z 3 = j. Using the fact that, 011 | e A ( A; ) t | < l 5 as wc ll as (38), 
we hnd 

< C(||too||c»+i+II«oIIc»+i; 


ipe 


—u^krt 
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Note the extra N + 1 powers of k come from the N + 1 antiderivatives in the Taylor 
Remainder formula. We need to estimate 


ife 


—iSTk 2 t 






where 

q + 2 1\ + 3/2 + 4/3 = N + 1 
We being by noting that, since k E h2 < , 


t 




Q+^l+^2+^3 


j^3(?+2/i +/ 2 +W+I 


p.3q+2l 


1+^2+W+l 


L 2 




^-0 r— +^2 H - ^3 


fc | g +3Zi + |; 2 +4i3 /^.xI'Z+'l + T 


^yl 9 + 2 Zi + f Z2+3Z3 ^ y 


< C 


+Z 2 +Z 3 


f 9+3/l +1 h +4^3 


f 9+2/l +1 12 +3^3 


where C is independent of za Therefore, since z/^ 


~ z/ 


-1 


tfe 




<7+/l+/2+Z3 


^, 39 + 2 / 1 +/2+AT+1 


L 2 


< c 

< c 

< c- 

l 

= c- 


^P^—b'T^t^.q+h+h+h 

+Z2+Z3 

( Ut t ) 


^|g+3/i + |Z2+4Z3 


^/f 9+2Zl + |Z2+3Z3 


§ <7+2 /1 +1 12 + 3/3 


1 _ 1 / 5 
4 2 V 2 


L 2 

9 + 3 ^ 1 + 1^2+4/3) 


^9+Zl+Z2+Z3 — ^ — § (f g+3/i + ^/2+4Z3^ 

,§ 9+2/1 + 1 Z2+3/3— ^ — ^ (I 9+3/1 + ^2+4/3) 

. N+l 1 
t 4 4 


K 
V 4 


where we used (42) in the last equality. 

Using a similar calculation, we can bound the L 2 norm of each j th derivative of this 
remainder term. One can show that for each integer triple l + s + r = j, we have 


_ N_ _ 1 
t 4 2 


‘9E«CH < C(||*„|| CJ + llf-ollc-,)—^75 


s+r 

2 


The proposition follows from the fact that s + r < j. □ 

Remark 3.6. The key point is that we can analyze the asymptotic behavior of w and 
v to any given order of accuracy 0{t ~ M ) (when t > O(^)) by choosing N (and hence 
m) sufficiently large and studying only the behavior of e~ UTk ^ t wjf )W and e~ UTk ^ t v{f w . 
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4 Decomposition of Solutions and Proof of the Main 
Result 


In this final section, we state and prove our main result. 

Theorem 4.1. Given any M > 0, let N > 4 M, and let m > N + 1/2. If the 
initial values wq,vo of © lie in the space L 2 (m), then there exists a constant C = 
C(m, N, wo, Do) and approximate solutions w app , v app , computable in terms of the 2N+ 
3 dimensional system of ODEs (14), such that 


H£,t) -ui apP (£,,T) || L 2 (m) + ||v(f,r) ~V app (f,T) || L 2( m ) < 


C 


-Mr 


N | m 
V 4^2 


for all t sufficiently large. The approximate solutions w app and v app satisfy equations 


(49) and (50) respectively. The functions 4>j(C) are the eigenfunctions of the operator 
C T (corresponding to diffusion with constant v T = v + ( in scaling variables) in the 
space Lf(m). The quantities a^r) and fkij) solve system (14) and have the following 


asymptotics, obtainable via a reduction to an N + 2-dimensional center manifold: 


KO)l < < 


C(Af,fc)e~t T 


C(N,k)e 


. fc+i. 


C(N,k)e~ 


_fc+2. 


C(N,k)e~ 


k+ 3. 


: k — 0 mod 4 
: k — 1 mod 4 
: k — 2 mod 4 
: k = 3 mod 4. 


(43) 


i aw i < < 


C(W,fc)e~^ T 


yk -\-1 


C(N,k)e~ 


M-i - 


/k -\-1 


C(N,k)e~ 


.k+2 n 


iyk -\-1 


C(N,k)e 


_k+ 3, 




: k — 0 mod 4 
: k — 1 mod 4 
: k — 2 mod 4 
: k — 3 mod 4. 


(44) 


Remark 4.2. As we will see in the course of the proof of the theorem, r > 0(log( ^ lo ^ t/ ^ )) 
(or equivalently t > 0(\ogv~ 1 )) will suffice for these estimates to hold. 


Proof of Theorem 4.1 We first concentrate on defining w app and v app and estab¬ 


lishing the error estimates in Theorem |4.1| this process will mainly use results from 
section 3. Recall the decomposition of w from section 3: 


w{k, t ) = ife UrkH (w% w + w\Z) + ^ge x ~ {k)t + w high . (45) 

The main results of section 3 essentially said w ~ f>e~ l " Tk ^ t w(f )W , with errors (measured 
in the ||| • ||| norm introduced in that section) either algebraically or exponentially 
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decaying. More precisely, using Propositions |3.3 and 3.4 we obtain 


| w — ipe 


-u T k 2 t-N 


W 


low III — 


< c 


JV I m N I 1 
V 4 2 t 4 ^2 


-L — u f 

H-—e s* 


z/ 


m+2 


where C is independent of u. For some t sufficiently large, we can “absorb” the 
exponentially decaying term into the algebraically decaying term; i.e. 


/m+2 


— »t 

e s < 


i 

JV I m JV.l • 
|/4 + 2f4 t 2 


We want to quantify how large t must be for the above inequality to hold. However, 
there are several other places in this section where terms of the form u~ p e~A t appear, 
which we wish to absorb into algebraically decaying errors. For this reason, we state 
and prove the following lemma: 

Lemma 4.3. Let A,M,E,p > 0 with v > 0 as before. Then there exists a constant 
C = C(M , A) > 0 such that for all t > log (u e ~ p ~ M ), we have the inequality 

1 1 

—e A < - 

up iAt M ' 

Remark 4.4. Note in particular that since r = log(l + t), the inequality t > 
£ log (u e ~ p ~ M ) essentially translates to r > 0(log( ^ lo +^ )). 


Proof of Lemma |4.3| We introduce a few new quantities to simplihy the notation: 
we let d = u p ~ e and we let a = v/A. Then the target estimate in the lemma reads 

t M e~ at < d. 

Now set f\(t) = t M e~ aM where 0 < A < 1 is fixed. Now, the target estimate in the 
lemma reads 

fx(t)e~ a ^- x)t < d. 

Using basic calculus, we find that the maximum value of f\ lies at t = and for 
t > g, we have f\(t) < (^) M ■ Therefore, if 


(- A f) e - a(1 - A)f < d, 

\a\e J 


we have the target estimate. The above inequality holds for 

M' 


t > 


log [d 


a( 1 — A) 

or, substituting a = u/A and d = u p ~ £ , we have 
A 


aXe 

~W 


t > 


1-A 


^ ^log(z/ £ p M ) + M ^log(Af) + log(^))) • 
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The time estimate in the lemma is just a less precise version of this inequality. This 


concludes the proof of lemma 4.3 


□ 


Next, we apply the lemma. Using the definition of w and inverting the Fourier 
Transform, we obtain, for t sufficiently large, 


Kx,t)-.F 1 ['ijj(k)e VTk2t w%J(x, 


< 


C 


N 

V 4 " 


-(1 + 1 ) 


-N/ 4 


(46) 


Proceeding, notice we can “drop” the cutoff function ^ in the above estimate with 
only an exponentially decaying penalty: due to the fact that 


\^{k)e 


—u T k 2 t-N 
w low 


— e 


-v T krt-N 


W, 


low I 


= \(m-l)e- rk2t 


w low I 


= 0 


for \k\ < which implies that 


I \8i ((m - l)e-^ t2 ‘wQ ||y < Te-P. 

From here on out, we will sometimes suppress the //-dependence of the constants for 
notational convenience. Proceeding, we define our approximate solution in x and t 
variables: 

jr-i \e- u T k2 tw% w \(x,t) = w app (x,t), 
which gives us the estimate 

III w(x,t) -w app (x,t) III < C(1 + t)~ N/4 . 


This is just estimate (|46|) without the cutoff function; it holds for t sufficiently large 

3, using scab 

w app (x,t) = 


as in Lemma T3 Therefore, using scaling variables and defining 

1 


\/l +1 


^appi^i '?’)) 


we have the estimate, which holds for r > O(l0g(M)), 


KM -w ap p(£,r) || L 2 (m) < 


C 


N | m 
Z/ 4 ‘ 2 


-K t 

e i . 


(This holds since the 111 • 111 and 11 • \ \L 2 (m) norms are equivalent in the way made precise 
at the beginning of section 3.) Using similar calculations, we have functions v app (x, t) 
and v app (£,r) satisfying 


and 


\\\v(x,t) ~ Vapp{x,t)\\ \ <C(l + t) N/4 

C _N 

v{t,T)-v app {£,T) || L 2 (m) < -^-^e T . 

V 4 + 2 
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This establishes the error estimates in 4T the remainder of the section is devoted to 
making more explicit the relationship between our approximate solutions w app , v app , 
and our center manifold calculations in £j2] 

Observe that 


'UJappyX 


,t) = [kie-^\x,t) 

n J’ 


3 = 0 
N 

E 

3=0 

N 


d jMM d3jr -i\ e -, T kH 


KM) 


y 9 3 k w(0,t) ^ 


3=0 


J KhF 


1 

g 4l Srjpt 


\J 47T Vxt 


Defining new scaling variables 


x 


and defining 


gives us 


f := r : = lo g(^) ; 


W a pp(x,t) ■ ^-Wapp^, t ), 


= E^De-^(0 oK ) 


i=o 

N 


n iy 




3=0 




(47) 


where the </>.,■(£) above are again the eigenfunctions of the operator on the space 
L 2 {m). 


We now show that the coefficients in (47) can be expressed in terms of the functions 


(afc(r), /3k(r)} from (j2j demonstrating that the leading order asymptotic behavior of 
the solution is determined by the center-manifold. 


First, recall from q3l formulas (34) and (36) that 


w(k,t ) = e UTk t w(k,t) + g{k)e 


Differentiating, we have 


a -t 


d J k w(k, t) = l {e UTk2t )diw(k,t)+dl(g(k)e x - t ) 

= + dl(g(k)e x ~ t ), 


1=0 
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where 1 is a polynomial of degree j — l. Setting k — 0, and substituting t = e r , 
we have 


d{w( 0 , e T ) 


E 

1=0 


C],ie 


^ f d l k w(0,e f ) + 0(e~^ f ), 


(48) 


where CJj = (j)i / T 2 P j *(0), and <9^(g(/c)e A “ t )|fc =0 is 0(e a eT ) since A_(0) = — v. We 
will proceed by computing the derivatives dj,w(0,t) in terms of the ay from (j 2 j 


Recall from £J2j formula (13), that we have the decomposition (using the orignal 
scaling variables £ and r) 


N 

w(£,r) = w c (£,r) + w s (£,r), w c (£, r) = aj( r )P?(0; w s = (w-w c ), 

3 =0 

N 

v(£, r) = u c (£, r) + u s (£, r), v c (£, r) = /3 j(t)^(£), = (v - v c ), 

3 =0 


and note the following. 

Lemma 4.5. J £ fc u; s (£, r)d£ = f £ fc u s (£, r)d£ = 0 for all k < N. 


Proof: We will prove the result for only, as the proof for v s is analogous. Note 
that 

N 

w s = w — P n w = w — Hj , w)4>ji 

3 =0 

and so ( H k ,w s ) = 0 for all k < N. We’ll proceed by induction on k. The k = 0 case 
follows because £° = 1 = Next, 


0 = ((H k+1 -H k ),w s ) = c k+1 / £ k+1 w s (f)df+^2c k / f k w s (f)df = c k+1 / £ k+1 w s (f)df 

J j = 0 J J 


by the inductive assumption. Since c k +1 7 ^ 0, the result follows. 
Using this lemma, we can compute 


□ 


d l k w(0,t) = 


d l k / e lkx w(x,t)dx 


I k=0 


4 / e ik ^w(l t )di 


fc =0 


(iVt+1) 1 J £ Z w(£,r)d£ 
(iVt + 1) 1 [ £V(£,p)d£ 


33 









for all l < N, and similarly for v. As a result, we have a relationship between d l k w(0, t) 
and the quantities a r : from i|2] 


N 


9 l k w( 0, t) = a r(log(l + *)) / X l (f> r ( 


X 


\/l + t 

v r=0 


y/l+t 


)dx, 


or equivalently 


N 


d l k w(0,e T ) = ^a r (log(l + e T ))(l + e r )5 f £ l MO d t 

r =0 d 


Inserting into (48), we obtain 


j N 

0,« f ) = J]q,e <i i i) ^ ar ( | og( 1 + e f ))(l + e , )U«Vrm + C l (e-5'') 

r =0 d 


1=0 


N 


N 


= e< 


^(Mi + qO^i + Osq, / + o(e-V). 


r=0 


1=0 


Therefore we can replace the coefficients in (47) and write 

N ( I N i _ , 

Wapp(i,f) = AAy Zl Q!? -(log(l + er ))Z]( 1 + e_r )' C '^ / £ 

j =0 V/ - '' ' r=0 1=0 J 


hdu 49) 


where C'J, ~ z/ J (see the line after (48)), and where we also have omitted an error 
term of 0{e~^ eT ) (This can be absorbed into the estimate in the original definition 
of w app by applying Lemma 4.3 with t = log(f).) Analogous calculations give us a 
similar result for v: 


N 


N 


J app 


(I, f) = £ ==] £ A-(log(l + e f )) £(1 + / t'Mm *(t).(50) 

3=0 ) r =0 1=0 ’’ J 


This completes the proof of Theorem 4.1 


Remark 4.6. Note that Theorem 4-1 is stated in terms of the original scaling vari¬ 
ables f and t. Since r = log(l + e T ), errors in f, for large f, are equivalent to errors 
in r, for large r. 


A Appendix: Convergence to the center manifold 


The purpose of this appendix is to show that the center manifold constructed in 
section 2 attracts all solutions. We know that we can construct the invariant manifold 
for the equations (21) (for a! k and b' k ) globally since we have explicit formulas which 


34 














hold for all values of a k and r/. In this note we show that any trajectory will converge 
toward the center manifold on a time scale of 0(1/u). 

Write bk = Bk+hk(ak- 2 ,..., do, rf). We’ll prove that for any choice of initial conditions 
Bk goes to zero like ~ e~ vt . 

First note that 

k—2 

b 'k = B ' k + Y ( d *t h k)a>i + (dv h k)r/ 

1=0,even 

k—2 ^ k—2 

= B'k - Y ( d a e hk)r)-di - Y (da e h k )vbe -2 - ( d v h k )rf 

1=2,even t=2,even 

k—2 k—2 ^ 

= B' k — rj Y (d ae h k )B f _ 2 + Y (^o^fcX-^a* - rjhe- 2 ) ~ (d v hk)r] 2 

t=2, even 1=2,even 


Thus, using the equation for b' k in (21) we have 

B k — V "Y; (b) ae hk)Bi _ 2 + (u + r]—)B k + -^-B k -2 

t=2,even 

1c 71 

= -(v + V7;) h k + -Wfc_ 2 
2 v A 

l 2 ij 

- Y ( d ai hk)(-r)-di - rjhe - 2 ) - ( d v h k )rf - ~jh k -2 

1=2,even 

The key observation is that the terms on the right hand side precisely represent the 
invariance equation that defines h k (and hence they all cancel), leaving the equation 

Jc 2 2 

B 'k = ~( u + r l^)Bk ~ ~B k -2 + T) Y ( d a e h k ) B e-2- ( 51 ) 

i=2,even 


We now see that the system of equations for B k is homogeneous, linear, upper- 
triangular (but non-autonomous), and hence can be analyzed inductively. We’ll show 
that 


B k (t) | < 


C(N) t 

u k/2 e 


for t > K We’ll only prove this for the even-indexed subsystem; the proof for the 
odd-subsystem is analogous 

k = 0 reads B' 0 = —vBq, which implies B 0 ~ e 
induction argument: assume for j — 0, 2,... £;, that 


Notice that the base case, k — 0, holds since (51) for 
~ ut Now let’s proceed with the 


BM < 


cm 

iyj/ 2 
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Next, we write the equation for B k +2 from (51) with a key difference in the way the 
last term is written: 


K+ 2 = - " B k +2 + V 


k + 2 


B 


k+2 


- —Bk+- c 

V 7 / ^' 


t= 1 


fc+2 ZLr ('co’i 

'fc+ 2-2i 


This last sum is obtained from reindexing £ and using the fact that the formula for 
hk in Proposition 2.6 implies that (d ae h k ) consists of a single term. Let’s proceed by 


noting that the equation 


y' = ~ {vy + arjy) 


has the exact solution 


y = e -^ t - to \l + t)- a {l + t 0 ) a 


To derive this solution, it may help to recall that y = Applying this to (52) with 


a = and using Duhamel’s formula, we obtain 


k+2 


-1 


Bk+ 2 (t) — e to \l+t) 2 (1 + t 0 ) 2 14 + 2 ( 0 ) + D^. +2 (t) + (53) 


where the Duhamel terms satisfy 


Dttut) 


c 


2£+l 


V 


e -i/(i-s)(i + + 5)^(1 + S ) 


\-l -1 


'to 


£=1 


(54) 

i/ 2 


Notice in the above Duhamel term, we have substituted, using the induction hypoth¬ 
esis \B k _ 2 f.{t)\ < fc -21 (we also assume t > t 0 > -). These Duhamel terms are 


the most slowly decaying terms in the solution formula (53). Proceeding, we simplify 


(54) and obtain (for all £), 




c 


V 2 

C 


+i+i 


e ut {l + t) ( fe 2 2 )^(i + t )2 £+1 — (1 + 4 ) 2 - 

1 1 (l + 4 )^ +1 ' 


i/t 


1/2 


+ 1 


i/(l + £)< zA(l + t 0 )* (1 + t)i +1 


Now since t > t 0 > -, we obtain 


C 


z/ 2 


and subsequently we obtain, for t > t 0 > -, 


C 


\B k+2 {t)\ < —j^e 


—vt 


V 2 


as desired. 


□ 
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